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Abstract. We study the planar and scalar reductions of the nonlinear Lindemann mechanism of unimolecular 
decay. First, we establish that the origin, a degenerate critical point, is globally asymptotically stable. Second, we 
prove there is a unique scalar solution (the slow manifold) between the horizontal and vertical isoclines. Third, we 
determine the concavity of all scalar solutions in the nonnegative quadrant. Fourth, we establish that each scalar 
solution is a centre manifold at the origin given by a Taylor series. Moreover, we develop the leading-order behaviour 
of all planar solutions as time tends to infinity. Finally, we determine the asymptotic behaviour of the slow manifold 
at infinity by showing that it is a unique centre manifold for a fixed point at infinity. 
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1. Introduction. A unimolecular reaction occurs when a single molecule undergoes a chemical 
change. For unimolecular decay (or isomerization) to occur, a certain amount of energy must be 
supplied externally, namely the activation energy. For some time, there was debate concerning just 
how the molecules became activated. Frederick Lindemann suggested [15] in 1922 that unimolecular 
decay involves two steps, namely the activation/deactivation by collision step and the reaction step. 
Cyril Norman Hinshelwood made further contributions |12) to the Lindemann model in 1926 and, 
consequently, the Lindemann mechanism is occasionally referred to as the Lindemann-Hinshelwood 
mechanism. For general references on unimolecular reactions and the Lindemann mechanism, see, 
for example, [H 1 [HI [15] . 

Suppose that the reactant A is to decay into the product P. Then, according to the (nonlinear, self- 
activation) Lindemann mechanism, A is activated by a collision with itself producing the activated 
complex B. This activation can also be reversed. The complex then decays into the product. 
Symbolically, 

A + A^A + B, b'^P, (1.1) 

k-i 



where ki, and k2 are the rate constants. 
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1.1. Differential Equations and Common Approximations. Using the Law of Mass 
Action, the concentrations of A and B in (jl.ll) satisfy the planar reduction 

da ^ db ^ 

— — k^iab — kia , — — = fcio — /c_iao — K20, (1-2) 
dr dr 

where r is time. The traditional initial conditions are a(0) = ag and 6(0) = 0. However, we will 
allow the initial condition for b to be arbitrary. Note that dp/dr = and (traditionally) p{0) = 0. 
Since the differential equations (jl.2l) do not depend on the differential equation for p, we need only 
consider the differential equations for a and b. 

There are two common approximations for the planar reduction. The Equilibrium Approximation 
(EA) and the Quasi-Steady-State Approximation (QSSA), which have proved successful for the 
Michaelis-Menten mechanism of an enzyme-substrate reaction jl7| . have also been applied to the 
Lindemann mechanism. See, for example, §2.2 of [TT] and pages 122-126 and 313-317 of [H]. 
These approximations are frequently employed to simplify more complicated networks in chemical 
kinetics which may involve, for example, inhibition or cooperativity effects. For the EA, one assumes 
da/dr « for sufficiently large time. This yields 

&(r)«^a(r). 

The QSSA, on the other hand, assumes db/dr w for sufficiently large time. This yields 

fcia(T)2 



6(t) 



k2 + fc_ia(r) 



It will be useful for us to convert the planar reduction to dimensionless form. Define 

t := k2T, X := \ — ] a, y := \ — ] b, and e :— - — > (J, 
V fc2 / V K2 / fci 

which are all dimensionless. Thus, t, x, and y are, respectively, a scaled time, reactant concentration, 
and complex concentration. Moreover, the parameter e > measures how slow the deactivation of 
the reactant is compared to the activation. Traditionally, one may want to consider s to be small. 
In our analysis, the size of e does not matter. 

It is easy to verify that, with the above rescaling, the planar reduction (|1.2p becomes 

X ~ —x^ + exy, y — x^ — {1 + ex) y, (1-3) 

where ' = d/dt. Observe that the system (|1.3p is a regular perturbation problem. Occasionally, we 
will need to refer to the vector field of this planar system. Hence, define 

where x := {x, y)'^ ■ Moreover, we will be working with the scalar reduction 
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where ' = d/dx, which describes solutions of the planar reduction (I1.3P in the xy-plane by suppress- 
ing the dependence on time. We will need to refer to the right-hand side of the scalar reduction. 
Hence, define 

fi^.y)--'^^^^. (1.6) 

—x^ + exy 



Remark 1.1. The function f(x,y) can be written f(x,y) = g2{x,y)/gi{x,y), where gi{x,y) and 
92{x,y) are the components of the function g(x) given in (|1.4p . Note the use of the row vector 
{x, y) in the arguments of gi and 52 as opposed to the column vector x. To alleviate notational 
headaches that arise from competing conventions involving row and column vectors, when there 
will be no confusion we will use the notation appropriate for the given situation. 



1.2. Discussion. The Lindemann mechanism has been explored mathematically by others. 
For example, the planar system (|1.3I) has been treated as a perturbation problem in [22l [24] . 
Furthermore, Simon Fraser has used the Lindemann mechanism [5J llOj as an example in his work 
on the dynamical systems approach to chemical kinetics. Finally, properties of the Lindemann 
mechanism have been explored mathematically in [7]. 

The focus of this paper is the detailed behaviour of solutions to the planar reduction (jl.Sp in phase 
space. That is, we perform a careful phase-plane analysis to reveal important details that a common 
phase-plane analysis would miss. Equivalently, we are studying solutions of the scalar reduction 
(|1.5p . It is worth reiterating that our analysis does not depend on the size of the parameter e 
(which is traditionally treated as being small). In ^ we present the basic phase portrait in the 
nonnegative quadrant. Moreover, we establish that the origin is a saddle node and is globally 
asymptotically stable with respect to the nonnegative quadrant. In ^31 we describe the isocline 
structure which we exploit in later sections. For example, the isocline structure plays an important 
role in determining the concavity and asymptotic behaviour of solutions. In 21 we prove that there 
is a unique slow manifold M between the horizontal and vertical isoclines. To this end, we use a 
nonstandard version of the Antifunnel Theorem. In ^J5l we determine the concavity of all solutions, 
excluding the slow manifold, in the nonnegative quadrant by analyzing an auxiliary function. In f}6l 
we use the Centre Manifold Theorem to show that all scalar solutions are given by a Taylor series 
at the origin. Moreover, we establish the leading-order behaviour of planar solutions as t ^ 00. 
This is nontrivial due to the fact that the origin is a degenerate critical point. In 53 we show that 
all planar solutions must enter and remain in the region bounded by the horizontal isocline and 
the isocline for the slope of the slow manifold at infinity. In SJSJ we single out properties of the 
slow manifold. These properties include concavity, monotonicity, and asymptotic behaviour at the 
origin and at infinity. Finally, in [JHl we state some open problems. 

2. Phase Portrait. A computer-generated phase portrait for the planar reduction p.3p . re- 
stricted to the physically relevant and positively invariant nonnegative quadrant 5, is given in 
Figure mU In this paper, we will develop precise mathematical properties of the phase portrait. 
Equivalently, we develop results on solutions of the scalar reduction (|1.5I) . 

The horizontal and vertical isoclines for the planar system (|1.3p . which are found by respectively 
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X ' = - + epsilon x y epsilon - 2.0 

y ' = - y - epsilon x y 




Fig. 2.1. A phase portrait for (11.31 1 for parameter value e = 2.0 along with the isoclines. 

setting y = and x = 0, are given by 

2 

y = H{x):^-^ and y^V{x):^-. (2.1) 

1 + ex £ 

The QSSA corresponds to the horizontal isocline (the quasi-steady-state manifold) and the EA cor- 
responds to the vertical isocline (the rapid equilibrium manifold). Observe that H{Q) = = V'(O), 
both H and V are strictly increasing, and V{x) > H{x) for all a; > 0. It appears from the phase 
portrait that the region between the isoclines, 

ro:={ (x,2/) : x>0, H{x) <y <V{x) } , 

acts like a trapping region for (time-dependent) solutions of the planar reduction. Moreover, the 
origin appears to be globally asymptotically stable. 

Theorem 2.1. Consider the planar system (jl.3p . 

(a) The region Tq is positively invariant. 

(h) Let x{t) he the solution with initial condition x{Q) = Xq, where Xq (z {x Cz S : x> 0}. Then, 

there is a t* > such that x{t) G Fq for all t > t* . 
(c) Let x(t) be the solution with x{0) = Xq, where Xq Cz S . Then, x{t) as t ^ oo. 
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Proof: 

(a) It follows from the definition (|1.4|) of the vector field g that g • < along V and H, where 
u is the outward unit normal vector. Thus, solutions cannot exit Fq through the horizontal 
or vertical isoclines. Furthermore, solutions cannot escape from Fq through the origin since 
solutions do not intersect. Hence, Fq is positively invariant. 

(b) Wc will break the proof into cases. 

Case 1: {xo,yo) E Fq. Since Fo is positively invariant, x(t) e Fq for all t > 0. 

Case 2: xq > and yo > V{xo). Suppose, on the contrary, that x(t) does not enter Fq. It 
follows that y{t) > V{x{t)) for all t >0. Using the differential equation (jl.3p . we know 
x{t) > and y{t) < for all t > Q. Now, we see from the definition (|1.6p of the function / 
that 

^ = fixit), y{t)) = -1 tt4^ TT^ < -1 for ah t > 0. 

•'vw:tfvy; ex{t)y{t) - x{ty 

Note that exy — x^ > since y > V{x) = x/e and x, y > 0. Thus, 

y{s) < —x{s) for all s > 0. 
Integrating with respect to s from to i and rearranging, we obtain 



y{t) < yo - [x{t) - xq] for all t > 0. 



Let {xi, V{xi)) be the point of intersection of the vertical isocline y = V{x) and the straight 
line y = yo — (x — xq). Obviously, xi > xq. Since x(t) is monotone increasing and bounded 
above by xi, we see that there is an a; G [a;o, Xi] such that x{t) — )■ a; as t — >■ oo. Similarly, 
since y{t) is monotone decreasing and bounded below by V{xo), we see that there is a 
y € [y(a;o),2/o] such that y{t) — )■ y as t — oo. Thus, the w-limit set is w(a;o,yo) = {(^il/)}- 
Since uj^XQ^yo) is invariant and (0,0) is the only equilibrium point of the system, a; = 
and y = 0. This is a contradiction. 
Case 3: xq > and < yo < H{xa). This case is proved in a manner similar to Case 2. 
(c) If Xq — 0, the solution of (|1.3p is x(t) = (0, yoe^*) . This clearly satisfies x{t) — ^ as t — > oo. 
Thus, we can assume a;o > and, by virtue of part (b), we can assume further that (xq, yo) G Fq. 
It follows from the differential equation (|1.3p and the fact that Fq is positively invariant that 
x{t) < and y{t) < for all t > 0. Since both x{t) and y{t) are decreasing and bounded below 
by zero, by the Monotone Convergence Theorem we know that there are x and y such that 
x{t) — X and y{t) y as t ^ oo. Thus, the w-limit set is lo{xo, yo) = {{x,y)}- Since uj{xo, yo) 
is invariant and (0, 0) is the only equilibrium point of the system, x = and y = 0. 

□ 



The Jacobian matrix at the origin for the planar system (|1.3p is diag(0,— 1). Thus, the origin is 
a nonhyperbolic fixed point. The Hartman-Grobman Theorem, unfortunately, cannot be applied 
here. Using Theorem 65 in §9.21 of [Ij, the origin is a saddle node which consists of two hyperbolic 
sectors and one parabolic sector. As we will effectively show later, S is contained in the parabolic 
sector. 
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Fig. 3.1. Graph of the function K{c). 



3. The Isocline Structure. The horizontal and vertical isoclines, along with all isoclines 
between them, will be very useful. If we solve /(a;, y) = c for y, we obtain y = F(x, c), where 



and 



j^(a:,c):^ , x^-e-'Kic), (3.1) 

K [c) + ex 

K{c):^-^, (3.2) 
1 + c 

That is, y = F{x, c) is the isocline for slope c. Figure 13.11 gives a sketch of K. Note that each 
isocline, for c € R\ {—1}, has a vertical asymptote at x = —e~^K{c). 

Remarks 3.1. 

(i) The interior of the region Fq corresponds to < c < oo and G < K{c) < \. 

(ii) Two exceptional isoclines are y = V{x) (the vertical isocline) and y = which correspond, 
respectively, to 

X 

lim F(x,c) = — and lim F{x,c)—Q. 

c— ^oo s c— > — 1 

Claim 3.2. Let c e IR\ { — 1} and let w{x) :— F{x, c) he the isocline for slope c. Then, the derivative 
of w satisfies 

lim w'{x) =e"^ (3.3) 

x— f oo 

Furthermore, w is concave up for all x > -~e~^K(c) and satisfies the differential equation 

x^w' + w{ew- 2x) = 0. (3.4) 
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y 




y = v{x) 



X 



increasing c 
(from — oo to —1) 

c switches from 

+00 to —00 



increasing c 
(from —1 to +00) 



Fig. 3.2. Sketch of the isocline structure of 111.51 1. The isoclines above the vertical isocline have zero slope along 
the line y = 2x/e. 



Proof: The proof is straight-forward and omitted. 



□ 



Remark 3.3. The vertical isocline satisfies the limit p.3p and the differential equation p.4p . The 
isocline w{x) — (the isocline for slope —1) also satisfies the differential equation but does not 
satisfy the limit. 



The isocline structure is sketched in Figure 13.21 We will often appeal to the isocline structure. 
For example, if a scalar solution of (jl.Sp is above the line y = and below the horizontal isocline 
y = H{x), we know that —1 < y'{x) < 0. 



4. Existence and Uniqueness of the Slow Manifold. It appears from the given phase 
portrait, Figure 12. 1[ that there exists a unique solution to (|1.5I) that lies entirely in the region 
between the horizontal and vertical isoclines. To prove this, we will need to use a nonstandard 
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version of the Antifunnel Theoreni. See, for example, Chapters 1 and 4 of [T5] . 

Definition 4.1. Let I = [a, 6) or I = (a.b) be an interval (where a < b < oo) and consider the 
first-order differential equation y' ~ f{x,y) over I. Let a,/3 G C^(/,]R) be functions satisfying 

a'{x) < f{x,a{x)) and f{x,l3{x))<l3'{x) for all x e I. (4.1) 

(a) The curves a and (5 satisfying (|4.ip are, respectively, a lower fence and an upper fence. If there 
is always a strict inequality in (j4.ip . the fences are strong. Otherwise, the fences are weak. 

(b) If /3{x) < a{x) on I, then the set T is called an antifunnel, where 

r:={ {x,y) : xel, f3{x) < y < a{x) }. 

Theorem 4.2 (Antifunnel Theorem, p. 196 of |13|). Let T be an antifunnel with strong lower and 
upper fences a and (3, respectively, for the differential equation y' = f[x,y) over the interval I, 
where I — [a, oo) or I = (a, oo). Suppose that there is a function r such that 

r{x)<^ix,y) for all {x,y)eT and lim "(^) " ^(^) _n 



dy ' ' x-^co exp(^f^ r{s) ds) 

Then, there exists a unique solution y{x) to the differential equation y' = f(x,y) which satisfies 
I3{x) < y{x) < a{x) for all x G /. 

Remark 4- 3- The standard version of the Antifunnel Theorem applies to antifunnels that are 
narrowing. That is, where a{x) — /3{x) — > as a; — oo. This version applies to, for example, the 
Michaelis-Menten mechanism [J]. 

4.1. Existence-Uniqueness Theorem. We want to show that there is a unique scalar so- 
lution that lies entirely in the region Tq. However, the vertical isocline is not a strong lower fence 
since "/(a;, V{x)) = oo." This turns out to be a fortunate obstacle. 

Suppose that we want the isocline 

x^ 

w{x) := , < r < 1 

r + ex 

to be a strong lower fence for the differential equation (|1.5p for all a:: > 0. Note that the condition 
on r restricts the isocline to being between the horizontal and vertical isoclines. Note also that 

f{x,w{x))=K^^{r) for all a; > 0. 

Since w is concave up and satisfies the limit p. 31) . we know that w'{x) < e^^ for all a; > 0. Hence, 

e^^<K^^{r) w'{x) < f{x,w{x)) for all a; > 0. 

Since we want the isocline that will give us the thinnest antifunnel, we choose 
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Note that a{x) is the isochne for slope e ^ and 



l + s 

Hence, define the region 

Ti { (x, y) : x> 0, H{x) <y< a{x) } . 

Theorem 4.4. 

(a) There exists a unique solution y = M{x) (the slow manifold) in Ti for the scalar differential 
equation p.5[) . 

(b) The solution y = Ai{x) is also the only solution that lies entirely in Fq. 
Proof: 

(a) We have aheady estabhshed that a is a strong lower fence. To show that if is a strong upper 
fence, observe 

f{x, H{x)) = < H'{x) for ah a; > 0. 

Moreover, a{x) > H{x) for all a; > 0. By definition, Fi is an antifunnel. 
Suppose that x > and H{x) < y < V{x). Using ([2Tl) . 

9 2 

£x , ^ ex 

< sy < X and < x ~ ey < x . 

I + ex 1 + ex 

Rearranging, we have 

„ 1 + ^2: 1 , „ /I \^ 1 
0< < and 0< -+e < ^. 

X x-ey \x J {x~ey) 



Thus, using the definition (|f .6p of / we have 

^{x, y) = ^2 > f- + > e"- 

Hence, we can apply the Antifunnel Theorem with r{x) := e^. To see why, consider that 



I r{x)dx = oo and lim [a{x) — H (x)] = ^r-. — 
Jo ^ (1 



+ e) 



and thus 



a{x) - H{x) 

lim , ^ , \ , \ = 0- 

K^oo expy^ r(s)ds) 

We can therefore conclude that there is a unique solution y = A4{x) to (jl.Sp that lies in Fi for 
aU a; > 0. 
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(b) Let 2/ be a solution in Fq lying below A4. Since Ai is the only solution contained in Fi, y must 
leave Fi through the horizontal isocline. 

Now, let y be a solution in Fg lying above A4. Suppose on the contrary that y never leaves Fg 
and thus M{x) < y{x) < V{x) for all x > 0. Consider the isocline w^x) := F (a;, 2e~^), which 
satisfies a{x) < w{x) < V{x). Since w is a strong lower fence, the proof of the first part of the 
theorem can be adapted to show that Ai is the only solution contained in the region between 
the isoclines H and w. Thus, there is an a > such that y{a) = w{a). Now, w'{x) < 
for all X > and y'{x) > 2e^^ if y{x) > w{x). It follows from a simple comparison argument 
that w{x) < y{x) < V{x) for all x > a. So, y{x) > y{a) + 2e^^ {x — a) for all x > a. This is a 
contradiction since y{x) > V{x) for sufficiently large x. 

□ 



Remarks 4-5- 

(i) We are referring to the unique solution between the horizontal and vertical isoclines as the 
slow manifold. However, all scalar solutions in Fg are technically slow manifolds (and, as it 
turns out, centre manifolds). This is because, as functions of time, the solutions approach the 
origin in the slow direction. 

(ii) There is no isocline w{x) such that w{x) > H{x) and w{x) is a strong upper fence for all x > 0. 
To see why this is the case, suppose w{x) :— F{x, c), where c > 0, satisfies ■w'{x) > f{x, w{x)) 
for all x > 0. This is impossible, since /(x, w(x)) = c for all x > and w'(x) as x — > 0+. 

Proposition 4.6. Let y be a solution to (ll.5|l lying inside Fi for x G (0,a), where a > 0. Then, 
we can extend y{x) and y'{x) to say y{0) = and y'{0) = 0. 



Proof: Observe that 



( X } 

Mm a{x) = and lim = 0. 

a;-i.O+ x^0+ X 



Since 

< t/(x) < a(x) for all x G (0, a) , 

the Squeeze Theorem establishes y{0) = 0. Now, 

y{x) - y{0) a{x) ^ /n ^ 

< < tor all X e (0, a) . 

X — X 

Thus, by the Squeeze Theorem again as well as the definition of (right) derivative, we can say 
y'(0) = 0. □ 



4.2. Nested Antifunnels. The region Fi is the thinnest antifunnel for x > with isoclines as 
boundaries. However, we can find thinner antifunnels than Fi which are valid for different intervals. 
For an isocline w{x) :— F{x,c), where < c < s~^, to be a strong lower fence on an interval, we 
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Fig. 4.1. Graph of the function g(c) for arbitrary e > 0. 



need w'lx) < f(x,w{x)). Solving the equation = f{x,w{x)), as we shall see, gives x = ^(c) 

where 



K{c) 



(4.3) 



Note that 1 — ec > 0. We will quickly establish a few properties of ^(c). See Figure HTTl for a sketch 
of the function. 

Claim 4.7. 

(a) The function £^{c) satisfies 



lim ^(c) — 0, lim ^(c) =00, ana 



C{c) > for all c G (0,e"^) . (4.4) 



(b) The function (^{c) is analytic for all c G (0,e Furthermore, (^{c) has analytic inverse ^ ^(x) 
defined for all x > 0. 



Proof: The proof is routine, tedious, and omitted. 



□ 
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Proposition 4.8. Letce (0,e ^) andw{x) := F{x,c). 
(a) The isocline w satisfies 

{< f{x,w{x)), if Q<x<i{c) 
w'{x) 1= f{x,w{x)), if a; = C(c) 
y> f{x,w{x)), if a; > C(c) 

(h) The slow manifold satisfies 

w{x) < M(x) < a{x) for all x > 

Proof: 

(a) Note that f{x, w{x)) — c for all x > 0. If we set w'{x) — c, we obtain 

e (1 - ec) x^ + 2K{c) (1 -ec)x- cK{cf = 0. 

This has two roots, one negative and one positive. The positive root is given by x = ^(c) with 
^(c) as in (|4.3I) . It is a routine matter to confirm that w^x) < f{x,w{x)) when < a; < ^(c) 
and that w'{x) > f{x,w{x)) when x > ^(c). 

(b) It follows from the Antifunnel Theorem. 

□ 

5. Concavity. In this section, we will establish the concavity of all scalar solutions, except 
for the slow manifold, in the nonnegative quadrant. The concavity of the slow manifold will be 
established later. These results will be obtained by using an auxiliary function. Moreover, we will 
construct a curve of inflection points which approximates the slow manifold. 

5.1. Establishing Concavity. Let y be a solution to (jl.Sp and consider the function / given 
in (II. 6p . If we differentiate y'{x) — f{x,y{x)) and apply the Chain Rule, we obtain 

y"{x) = p{x, y{x))h{x, y{x)), (5.1) 

where 

P{x,y) := „ , ^ "2 and h{x,y) := x'^ f{x,y) + y {ey ~ 2x) . (5.2) 

x"^ (ey — X) 

The function p{x, y) is positive everywhere except along the vertical isocline and for x = 0, where 
it is undefined. We will be considering the functions h(x, y) and p{x, y) along a given solution y{x) 
so we will abuse notation by writing h{x) := h{x,y{x)) and p(a;) := p{x,y{x)). 

For a given x > with y{x) ^ V{x), it follows from (|5.ip and the fact that p{x) > that the sign 
of h{x) is the same as the sign of y"{x). Furthermore, if we differentiate h(x) with respect to x and 
apply (|5.ip . we see that the function h has derivative 



h'{x) = x^p{x)h{x) + 2y{x) [ey [x) - 1] 



(5.3) 
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Region 


Concavity of Solutions 


0<y <H 


concave down 


H < y < M 


concave up, then inflection point, then concave down 


M < y <V 


concave up 


y>v 


concave up, then inflection point, then concave down 



Table 5.1 



A summary of the concavity of solutions of 111.511 in the nonnegative quadrant. 



Remark 5.1. The function h cannot teU us anything about the concavity of solutions at a; = 0, not 
even by taking a limit. 

Claiivi 5.2. Let y he a solution to (jl.Sp and let > with y{xo) ^ V{xo). Consider the isocline 
through the point {xo,y{xo)), which is given by w{x) := F{x,y' (xq)) . Then, 

h{xo) = xo [y'{xo) - w'{xo)] . 

Furthermore, 

y"ixo) > <^=4> y'ixo) > w'{xo) and y"{xn) < y'ixo) < w'{xo). 



Proof: The first part follows from (I3.4p and (|5.2p . The second part follows from the first. □ 

The concavity of all solutions in all regions of the nonnegative quadrant can be deduced using the 
auxiliary function h and the following easy-to- verify lemma. Table 15.11 summarizes what we will 
develop in this section. 

Lemma 5.3. Let 1 be one of the intervals [a, b], (a, h), [a, b), and (a, h\. Suppose that cj) £ C(/) is a 
function having at least one zero in L . 

(a) If 1 — (a, b] or I = [a,b], then the function (f> has a right-most zero in I. Likewise, if I = [a,b) 
or I — [a, b], then the function (p has a left-most zero in I . 

(b) If 4> £ C^{I) and 4>' {x) > for every zero of (j) in I, then (j) has exactly one zero in I. 

Proposition 5.4. Let y be a solution to (|1.5p lying below H with domain [a,b], where < a < b, 
y{a) = 11(a), and y{b) = 0. Then, y is concave down on [a, 6]. 

Proof: We know y'(a) = and y'{b) = —1. Also, y{x) > and y'{x) < for all x e {a,b). Note 
that y < V{x) which implies ey — 2x < Q for all x E [a, b]. Let h be as in (|5.2p defined with respect 
to the solution y. Observe that 

h{a) = H{a) [eH{a) - 2a] < and h{b) = -6^ < 0. 

Observe also that 

h{x) — x^y'lx) + y{x) [£^(2:) — 2a;] < for all x £ (a, b). 
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Therefore, y is concave down on [a, &]. □ 

Proposition 5.5. Let y he a solution to (|1.5p lying above H and below Ai with domain (0, a], 
where a > and y{a) = H{a). Then, there is a unique Xi e (0, a) such that y" [xi) — 0. Moreover, 
y is concave up on (0,a;i) and concave down on (xi,cl\. 

Proof: Let h be as in (|5.2p defined with respect to the solution y. Now, we know y'{a) — and, by 
Proposition 14.61 we can extend y'{x) continuously and write 2/'(0) = 0. By Rohe's Theorem, there 
is an xi £ (0, a) such that y"(xi) = and hence h(xi) = 0. To show the uniqueness of xi, suppose 
that X2 G (0, a) is such that h{x2) = 0. Now, since H{x2) < 2/(2:2) < a{x2), by virtue of the isochne 
structure < y'{x2) < s~^. Moreover, we see from (|5.3|) that h'{x2) < 0. By Lemma ESI we can 
conclude X2 = xi. Finahy, by continuity we can conclude that h(x) > on (0, xi) and h{x) < on 
{xi,a] since h{a) — y{a) [ey{a) — 2a] < 0. □ 

Proposition 5.6. Let y he a solution to (|1.5p strictly between Ai and V with domain (0,a), where 
a > and y{a—) — V{a). Then, y is concave up on (0, a). 

Proof: Let h be as in (|5.2p defined with respect to the solution y. Define the sets 

/ := { X e (0, a) : M{x) < y{x) < a{x) } 

and 

J := { X e {0,a) : a{x) < y{x) < V{x) } . 

Since since y{a—) — V{a), we know that J 7^ 0. We will show separately that h{x) > for each 
X £ I and for each x £ J. 

Suppose that x G J. By virtue of the isocline structure, y'{x) > £~^. However, any isocline w 
satisfies w'{x) < e^^. Appealing to Claim [Ol h{x) > 0. 

If / = then we are done. Suppose that / 7^ and let X3 G J be fixed. We have already shown 
that h{x3) > 0. We need to show that h{x) > for all x € /. By continuity, it suffices to show 
that h has no zeros in I. Suppose, on the contrary, that this is not the case. By Lemma 15. 3[ h 
has a right-most zero xi e / C (0,X3]. Observe that y'{xi) < e~^, which follows from the fact that 
y(xi) < a(xi). Since h{xi) = and y{xi) > 0, (|5.3I) informs us h'{xi) < 0. Consequently, there 
is an X2 G (xi,X3) such that h{x2) < 0. Since h{x3) > 0, the Intermediate Value Theorem implies 
that h has a zero in (x2,X3), which is a contradiction since xi is the right-most zero. □ 

Proposition 5.7. Let y be a solution to (II. 5p lying above V with domain (0, a), where a> 0, 
y{0+) = 00, andy{a—) — V{a). Then, there is a unique Xi G (0,a) such that y" (xi) = 0. Moreover, 
y is concave up on (0,xi) and concave down on (xi,a). 

Proof: We know that 



Mm y'{x) — —00 and lim y'{x) = —00. 

x— >0+ a;— 
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Fig. 5.1. The two thick curves are curves along which solutions of 1)1. 5p have inflection points, for parameter 
value e = 0.5. The thin curves are the horizontal, a, and vertical isoclines. We will show later that the slow manifold 
lies between the lower thick curve and the middle thin curve. 

By continuity, there are bi, &2 € M such that < 6i < 62 < « and y'{bi) = y'(&2)- Thus by Rohe's 
Theorem, there is an xi G (0,a) such that y"{xi) — 0. 

To prove uniqueness, suppose that X2 G (0, a) is such that h{x2) = 0, where h is as in (|5.2|) defined 
with respect to the solution y. Now, we know from the isochne structure that y'{x2) < s^^. Using 
(j5.3l) . h'{x2) < 0. Since X2 was an arbitrary zero of /i, we can conclude using Lemma 15.31 that 
X2 = xi. Furthermore, since h'{xi) < 0, we can say that y is concave up on (0,xi) and concave 
down on (xi, a). □ 

5.2. Curve of Inflection Points. We know from Table 15.11 that solutions to the scalar 
differential equation (|1.5p can only have inflection points between H and M or above V. We can 
construct a curve of infiection points, between H and M, which is close to the slow manifold. 

It is easily verified that 



where h is as in ()5.2|) . Thus, there are three curves along which solutions have zero second derivative, 
given implicitly by 



One curve lies below the x-axis and is discarded. The other two curves, as expected, are in the 
positive quadrant. See Figure ISTTl 



h{x,y) 



e^y^ — (Sex) + (2a;^ — ex"^ — x) y + x"^ 



ey — X 



£^y^ - (Bex) + (2x^ - ex'^ - x) y + x^ = 0. 
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Recall that, for a fixed c G (O, e~^) , the isocline w{x) := F{x, c) switches from being a strong lower 
fence to being a strong upper fence at a; = ^(c) and y = F(^(c), c), where F is defined in (|3.ip and 
^ is defined in (|4.3|) . As it turns out, 



(5.4) 



will be a curve of inflection points between H and A4. Note that < F{x,c) < a{x) for all 

X > and c € (0,e~^), which follows from the isocline structure. Moreover, note < £,^^{x) < 
for all a; > 0. Thus, H{x) < y{x) < a{x) for all a; > 0. 

Claim 5.8. Suppose that xq > and H{xq) < yo < a(xo). Define the slope c :~ /(xo,2/o) '^i^d, 
isocline w{x) :— F{x,c). Then, the isocline w satisfies 



w'{xo) < 



> fixo,yo), if H{xo) < yo < yixo) 

^ f{xo,yo), if yo^y{xo) 

< f{xo,yo), if y{xo) < yo < a{xQ) 



Proof: Note that < c < e^^ and j/o = ^^(a^o)- We will only show the third case since the other 
two cases are similar. Assume that y^xo) < yo < a{xo)- Appealing to the isocline structure, we 
know fy{x,y) > if a; > and H{x) <y < a{x). Consequently, f{xo,yo) > fixo,yixo))- Since 
c = /(xq, j/o) and (,^^{xo) = /(xo, 3^(xo)), we can conclude c > ^^^(xq). Since ^ is strictly increas- 
ing, xq < ^(c). By virtue of Proposition 14.81 we can conclude w'{xo) < f{xo,yo)- □ 

Claim 5.9. The curve y — y{x) is analytic for all x > 0. 

Proof: We know that ^~^(x) is analytic and < ^"""^(x) < for all x > 0. Since F{x,c) is 

analytic if x > and < c < e^-^, we see from the definition (j5.4p that y{x) is analytic for all 

X > 0. □ 

Proposition 5.10. The function h, defined in (|5.2I) . satisfies 

r < 0, */ X > 0, Hix) <y< y{x) 
h{x,y)l^O, if x>0, y^y{x) 

[ > 0, if X > 0, y{x) <y < a{x) 



Proof: Let xq > and H{xo) < yo < a{xo) be fixed. Consider the slope c := /(xq, yo) and isocline 
w{x) :— F{x,c). We know from Claim that 

H^o, yo) = xo [/(xo, yo) - w'{xo)] ■ 
The result follows from Claim \5M □ 



Proposition 5.11. The curve y = y{x) satisfies 

H{x) <y{x) < M{x) for all x > 0. 
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Proof: We know already that H{x) < y{x) < a{x) for all x > 0. We know from our results on 
concavity (see Table ISTTj) that h{x,y) > if x > and A4{x) < y < a{x), where h is the function 
defined in (j5.2[) . By continuity, we can conclude h{x, A4{x)) > for all x > 0. It follows from 
Proposition 15.101 that y{x) < M{x) for aU x > 0. 

To establish a strict inequality, let h be defined along the solution y — A4{x). Assume, on the 
contrary, that there is an > such that h{xo) = 0. Using (|5.3I) . h'{xo) < 0. This contradicts the 
fact that h{x) > for aU a; > 0. □ 

5.3. Slow Tangent Manifold. The curve y = y{x) can be referred to as a slow tangent 
manifold (or as an intrinsic low-dimensional manifold) since it consists of the points for which 
the tangent vector for the planar system (|1.3p points in the slow direction. See, for example, 
[HI [TCI [T51 [5DJ [211 [22] . To see why the curve of inflection points and the slow tangent manifold are 
equivalent, first consider the general planar system x = g(x), where '= d/dt and g e C^(R^,R^), 
along with the corresponding scalar system y' — g2{x,y)/gi{x,y), where ' = d/dx. Consider the 
linearization matrix A(x) {ffy (x)}^ ^-^j^, where gij(x) :— dgi(jx.)/dxj, which has characteristic 
equation 

A^-tA + A = 0, where t:= 511+522 and A 511522 - 512521- 
For notational brevity, we are suppressing the dependence on x. 

Claim 5.12. Suppose 4A < t'^ and 512 7^ 0. Then, A has real distinct eigenvalues and associated 
distinct eigenvectors given, respectively, by 



t±Vt^-4A fl\ A± - 511 

A± := and v± := \ , where a± := . 

2 \<^±J 512 



Proof: The proof is routine. □ 

Proposition 5.13. Suppose that 51 7^ 0, 512 ^ 0, and 4A < at some fixed point (a, 6) and let 

y{x) be the scalar solution through {a,b). Then, y"{a) — if and only if g \\ 1)+ or g \\ u_ at {a,b). 

Proof: First, note that a± is the slope of the eigenvector v±. If we differentiate the scalar 
differential equation and manipulate the resulting expression, we obtain 

// ^ -512 [(52/51) - <^+] [(52/51) - g-] 
The conclusion follows. □ 



For the specific planar and scalar systems (II. 3p and (|1.5p . we have 

— 2.T + ey ex 
2x — ey —ex — 1 



A=[ ^""^^y T= - ie + 2)x + ey-l, and A^2x~ey. 
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To apply Proposition [2331 we need to verify that gi ^ 0, dgi/dy ^ 0, and > 4A in the relevant 
regions. Trivially, gi 7^ (except along the vertical isocline) and 512 > for a; > 0. To show that 
> 4A for a; > 0, observe 



- 4A = £2 



y : 



4£a; > 4£a; > 0. 



This establishes that 3^ is a tangent manifold. To establish that y is indeed a slow tangent manifold, 
we note that (as can be shown) A_ < A+ < 0, cr_ < < (t+, and > for every x g Fi. We 

have thus demonstrated the following. 

Proposition 5.14. Consider the planar system (jl.3p and the scalar system (11.51) . The curve of 
inflection points y — y{x) is a slow tangent manifold. 

6. Behaviour of Solutions Near the Origin. In this section, we establish the full asymp- 
totic behaviour of scalar solutions y{x) as a; — > 0+. Moreover, we will obtain the leading-order 
behaviour of planar solutions x(t) as t — > 00. 

6.1. Scalar Solutions. We will begin by attempting to find a Taylor series solution. Consider 
the differential equation (|1.5p . which can be rewritten 

sxyy' — x^y' ~ x"^ + y + exy = 0. (6.1) 

Assume that y{x) is a solution in Pq of the form 



y{x) = J2 bnx^ (6-2) 



Tl = 

for undetermined coefficients {b„}'^^Q. If we substitute the series (|6.2I) into (|6.ip and then solve for 
the coefficients, we obtain 

bo = 0, 61 = 0, &2 - 1, 63 = 2 - £, 

ri-2 

and 6„ = (n — 1 — £) — £ (n — m) bmbn-m for n > 4. (6.3) 

m=2 



We will use centre manifold theory to show that the series (|6.2p is fully correct for each solution 
inside the trapping region Pq. However, we must first show that each solution is a centre manifold. 
That is, we must show that each solution y{x) satisfies 7/(0) — and 2/'(0) — 0. Proposition 14.61 
already established that this is true for y{x) inside Pi. 

Proposition 6.1. Let y he a solution to (jl.Sp lying inside Pq for x G (0, a), where a > 0. Then, 
we can extend y{x) and y'{x) to say y{0) = and y'(0) = 0. 

Proof: These limits have already been established if y is the slow manifold or if y lies below the 
slow manifold A4. Hence, we will assume that 

M{x) < y{x) <V{x) for all xe{0,a). 
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Let c y' (a/2). We know from Table [01 that y is concave up on (0, a/2). Since M.{x) > for 
X e (0, a/2), we thus have 

< y{x) < F{x, c) for all x e (O, ^ a) , 
where F is the function given in (|3.ip . Note that 

F(x c) 

lim F(x,c)=0 and lim ^ ' ^ = 0. 

x~¥0+ x^0+ X 

It follows from the Squeeze Theorem that we can take y{0) = 0. Now, observe that 

^ ^ yjx) - yjO) ^ F{x,c) 
X — Q X 

Again by the Squeeze Theorem, we see that we can take y'(0) = 0. □ 

Theorem 6.2. Let y{x) he a scalar solution to (jl.Sp lying inside Fq and consider the coefficients 
{6„}^2 d'tven in (|6.3p . Then, 

oo 

y{x) ^ 6„x" as x 0+. 

n=2 



Proof: The Centre Manifold Theorem guarantees that there is a solution u{x) to ()1.5|) such that 

oo 

m(x) ^ ^ 6„a;" as a; — > 0+. 

Note that the Taylor coefficients of the series for u{x) must be {&n}J^2 since they are generated 
uniquely by the differential equation. Since y{x) is a centre manifold, it follows from centre manifold 
theory that 

y{x) - u{x) ^ 0{x'') as a; -> 0+ 

for any k € {2,3,...}. See, for example. Theorem 1 on page 16, Theorem 3 on page 25, and 
properties (1) and (2) on page 28 of [Q- The conclusion of the theorem follows. □ 

Remark 6.3. For analytic systems of ordinary differential equations for which the Centre Manifold 
Theorem applies, if the Taylor series for a centre manifold has a nonzero radius of convergence, then 
the centre manifold is unique. Since all solutions y to (jl.5l) lying inside Fq are centre manifolds, we 
can conclude that the Taylor series 'Yl^=2 ^nx" has radius of convergence zero. 

6.2. Planar Solutions. We can use the isoclines to extract the leading-order behaviour of 
planar solutions as time tends to infinity. 

Proposition 6.4. Let x{t) be the planar solution to (|1.3p with initial condition x(0) = Xq, where 
Xo e {xe S : a; > 0}. Then, 
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Proof: Let c > be fixed and arbitrary. We know from Theorem 12. 1[ Table 15.11 and the isochne 
structure that there exists a T > such that 

Hix{t)) < y{t) < F{x(t),c) for all t > T, (6.4) 

where F is given in (|3.ip . Using (II. 3|) . (|2.ip . (|3.ip . and (|6.4p . we can see that x{t) satisfies 



x{tf 



1 + £x{t) 



< xit) < for all t > T, 

- ^ ^ - 1 + bx{t) - ' 



(6.5) 



where b :— e/K{c) and K is the function defined in p.2p . Note that the solution of the initial value 
problem 



1 + au 

where a, mq > and to ^ are constants, is 
u(t) = (p{t;a,to,Uo) 



— exp^ 



exp 



■t-tp- 



where W is the Lambert W function f6^ . A simple comparison argument applied to ()6.5p establishes 
ip{t;e,T,x{T)) < x{t) < Lp{t;b,T,x{T)) for all t>T. (6.6) 



A standard property of the Lambert W function is 

W{t) = \n{t)-\n{\n{t))+o{\n{\n{t))) as t ^ oo. 
Consequently, it can be shown 



W{e*) ^t-\n{t)+o{\n{t)) and 



1 



T^(e*) t t2 
With a little manipulation, it can be verified that 



1 ln(t) f\n(t) 



X 1 Ht) fHt)\ 

ip[t;a,to,uo) = J + a—^ + o[—^ j as t 



■ oo 



for any a, mq > and ^ 0. It follows from (16. 6p that 



lim inf 



x(t) - - 



4-1 1 



ln(i)J 



> e and lim sup 



x{t) 



as t — ?► oo. 



In(<) 



< b. 



Since c was arbitrary with K(c) — 1 and & — > e as c — > 0^, 



lim 



x(t) 



ln(i) 
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This yields the desired conclusion for x{t). The conclusion for y{t) follows from Theorem 16.21 □ 



It is possible to derive the expression for x{t) in Proposition 16.41 without appealing to the isocline 
structure and concavity. To achieve this, we will note that x{t) satisfies the integral equation 



2:(s) 



ds 



(6.7) 



and then twice utilize the following easy-to-verify lemma. 



Lemma 6.5. Let a G M be a constant and let /, g : [a, oo) — > M fee nonnegative, integrable functions 
such that f{t) = g{t) + o{g(t)) as t ^ oo. If G{t) is an antiderivative of g{t) such that G{t) — >■ oo 
as t — > oo, then f{s) ds = G{t) + o{G{t)) as t oo. 



We know from Theorems 12.11 and 16.21 that 



yit) 



x{t) — o(l) , y(t) = o(l) , and + 1 = 1 + o(l) as t ^ oo. 

xit) 



It follows from Lemma 16.51 that 



^0 



x{s) 



ds ^t + o{t) 



x{s) 



ds — o{t) as i — )■ oo. 



By virtue of the integral equation (|6.7I) , 



1 



x{t) 



i [1 + 0(1)] =^ x{t) = - + 



as i — oo. 



as t — > oo, 



To take this one step further, observe now that 

yW=l + o( 

x{t) t \t^ 
which follows from Theorem 16.21 and so by Lemma 16.51 we have 

/ ^74rfs = ln(i) + o(ln(<)) as t 
Jo ^is) 

By virtue of the integral equation (|6.7|) once again. 



1 



= t 



t 



t 



u{t) = -+e 



1 \n{t) f\n{t) 



t2 



as t — > oo. 



7. All Solutions Must Enter the Antifunnel. Earlier, in Theorem 12.11 we showed that 
all solutions x(t) to the planar system (jl.3|) . except for the trivial solutions, eventually enter the 
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trapping region Tq. Here, we show that Fi is itself a trapping region. 

Theorem 7.1. Let x{t) he the solution to (|1.3p with a;(0) = Xq, where Xq E {x<Z S : x > Q}. 

(a) There is a t* > such that x{t) G Fi for all t > t* . 

(b) Define the region 

F2 {{x,y) : x> 0, y{x) <ij< aix) } . 
Then, there is a t* > such that x(t) G F2 for all t > t* . 

Proof: 

(a) We know from Theorem 12.11 that x{t) eventuaUy enters and stays in Fq. Let y{x) be the 
corresponding scalar solution to (|1.5p . Then, we can say y'(0) = 0. Appealing to the isocline 
structure, this means that x(i) has entered Fi. Furthermore, since g • i^* < along the horizontal 
and a isoclines which form the boundaries of the region in question, we see that Fi is positively 
invariant. 

(b) It follows from Table [STTl Proposition l5.10[ and the previous part of the theorem. 

□ 

8. Properties of the Slow Manifold. In this section, we will highlight some properties of 
the slow manifold. 

Proposition 8.1. The slow manifold y = M{x) satisfies, for all x > 0, 

< H{x) < y{x) < M{x) < a{x). 

Furthermore, 

lim M{x) = 0. 

x-)-0+ 



Proof: The first part follows from Theorem 14.41 and Proposition 15.111 The second part follows 
from the Squeeze Theorem. □ 

Proposition 8.2. The slow manifold y = A4{x) is concave up for all x > 0. 

Proof: It follows from Propositions 15. lOl and 18. II that h{x, Ai{x)) > for all a; > 0, where h is the 
function defined in (15.21) . Since sgn(A^"(a::)) — sgn(h(x,Ai(x))), it must be that the slow manifold 
is concave up for all a; > 0. □ 



Proposition 8.3. The slope of the slow manifold y ~ M{x) satisfies 

< M'{x) < e^^ for all x > 0. 
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Furthermore, 

Urn M'{x)=G and lim M'{x)=e~^. 

X-S-0+ x->oo 



Proof: The first part is a consequence of Proposition 18 . 1 1 and tlie isocline structure. Tfie first limit 
is a special case of Proposition 16.11 To prove the second limit, let cG (0,e~^). It follows from 
Proposition 14.81 and the isocline structure that 

c<M'{x)<e^^ for aU x>^{c), 

where ^ is the function defined in (14.31) . Applying (|4.4I) and the Squeeze Theorem gives the second 
limit. □ 

Remark 8.4- The justification which Fraser provides in [9| (just before Theorem 1) that M.'{x) — >■ 
as a; — cx) is incorrect. The error is that the distance between the horizontal and vertical isoclines 
does not tend to zero as x tends to infinity. Thus, the asymptotic behaviour of M' [x) need not be 
the same as the asymptotic behaviour of H'{x) and V'{x). 

Proposition 8.5. Asymptotically, the slow manifold can be written 

oo 

M{x) "^hnx" as a;^-0+, 

n=2 

where the coefficients {&n}^2 '^''^ '^'^ (|6.3p . 

Proof: Since the slow manifold is contained entirely in Pq, we can apply Theorem 16.21 □ 

Corollary 8.6. The slow manifold satisfies 

M{x) ^ H{x) +0{x^) as x ^ 0+. 
Moreover, this statement would not be true if we replace H{x) with any other isocline F{x,c). 

Proof: It follows from a comparison of the asymptotic expansions for A4{x), H{x), and F(x, c). □ 

Now we will establish the full asymptotic behaviour of the slow manifold at infinity. First, we will 
extract as much information as possible from the isoclines. Second, we will attempt to find a series 
in integer powers of x. Third, we will prove definitively that the resulting series is indeed fully 
correct. 

Let c e (0,e^^). We know from Proposition 14.81 that 



F{x,c) < M{x) < a{x) for ah x > ^(c). 
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where F is defined in p.ip and ^ is defined in (I4.3p . Note tliat 



F a:,c = V^ + O - and a a; = 



\ — \ as a; — > oo . 



Since 



F(a;,c) - - < M[x) - - < a{x) - - for all x > ^(c), 

e £6 



we can conclude 



lim inf 



M{x)-- 



> — and limsup 



M{x) - - 

£ . 



< 



£(! + £)■ 



Since c e (0,£ ^) is arbitrary and 



lim 



K{c) 



we have 



lim 

a:— J-oo 



X(a;) 



c^(e-i)- £2 £(! + £)' 



and A^(x' ^ 



£(! + £) 



£ £(! + £) 



o(l) as a; — )■ oo. 



Assume that we can write 



M{x) ^ p„a; " 



for undetermined coefRcients {Pn}°^^_i- Of course, we expect p_i = £ ^ and po = — e ^ (1 + £) ^ . 
Now, write the differential equation (|1.5p as 

£a;?/?/' — x^y' — a;^ + y + £xy = 0. (8-2) 

If we substitute (18.11) into (18.21) and solve for the coefficients, we obtain 



P-i — Po — 77—. — 7, 

£ £(! + £) 



and pn = - 



l + £ 



_i - £ ^ (n - m) p„i-iPn- 



m—1 



for n > 1. 



(8.3) 



Proposition 8.7. Asymptotically, the slow manifold can be written 

oo 

M{x)^ Pux"" 

n=-\ 

where the coefficients {pnS^=^\ '""6 



as a: -> oo, 
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Proof: To prove the result, we will apply the Centre Manifold Theorem to a fixed point at infinity. 
Consider the change of variables 

X x^^ and Y y ~ r{x), where r{x) :— p-ix + po + Pix^^, 

with the coefficients p^i, pa, and pi being given in ()8.3|) . Differentiate the new variables with 
respect to time and use the differential equation ()f .3|) to obtain the system 

X = -X^gi {X-\r{X-^) +Y) , 

Y = -r' (X-i) 51 {X-\r {X-^) + F) + {X-\r (X"!) + Y) , 

where gi and (72 are as in ^1A\ . This system is not polynomial but there is no harm in considering 
the system 

X^-XV {X-\r{X-^)+Y) , 

Y = X [-/ (X-i) 51 {X-\r (X-i) + r) + 52 {X~\r {X-') + Y)] , (8.4) 

which is polynomial. This is because the resulting scalar differential equation will be the same. The 
system at hand, while messy, is in the canonical form for the Centre Manifold Theorem. Note that 
the eigenvalues of the matrix for the linear part of this system are and — (f + e). We know from 
centre manifold theory that there is a C°° centre manifold Y = C{X) which, we claim, must be the 
slow manifold. 

For the scalar differential equation in the original coordinates, all other solutions except for the slow 
manifold leave the antifunnel Fi. To establish that the slow manifold in the original coordinates is 
the same as the centre manifold in the new coordinates, we need only show that Y = C{X) is the 
only scalar solution in the new coordinates which is o(l) as X ^ . 

Observe that the Y-axis, is invariant. Moreover, the fixed point (X, Y) = (0, 0) is a saddle node (or 
a degenerate saddle). The physically relevant portion of the phase portrait, namely X > 0, consists 
of two hyperbolic sectors, one with the positive K-axis and the centre manifold as boundaries and 
the other with the negative F-axis and the centre manifold as boundaries. See Figure 18.11 This 
can be shown using techniques in §9.21 of [1 (in particular Theorem 65 on page 340). Therefore, 
C{X) is the only scalar solution in the new coordinates which is o(l) as X — > 0+. It follows that 
the centre manifold is indeed the slow manifold. 

By the Centre Manifold Theorem, in the new coordinates the slow manifold can be written 

00 

X(X)-^p„X" as X^0+, 

for some coefficients {Pn}^2- Upon reverting back to original coordinates and observing that the 
coefficients in (|8.3p are generated uniquely from the differential equation, the conclusion follows. □ 

Corollary 8.8. The slow manifold satisfies 

M{x) — a{x) ^ o(^^i^ as x ^ 00. 
Moreover, this statement would not he true if we replace a[x) with any other isocline F{x,c). 
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X ' = 1/2 x^- 1/4 yi'-ify 

y' = -2y-3/8x^ + 1/16x^+1/4 x^ y - 




Fig. 8.1. A phase portrait for iSAt for e = 1.0, 



Proof: It follows from a comparison of the asymptotic expansions for M{x)^ a{x)^ and F{x^ c). □ 



9. Open Questions. It would be nice to extend Proposition 16.41 to include more terms. In 
particular, it is desirable to have the lowest-order term which depends on the initial condition. For 
x(t), it is expected that the initial condition xq first appears in the l/t^ term since this is the case 
when e = 0, which has 



x{t) 



Xq 



1 + Xot 



E 

n=l 



(-1) 



n+1 



as t ^ oo. 



More generally, we would like to develop an iterative procedure to extract as many terms as possible 
from the asymptotic expansion of a solution of a nonlinear differential equation which approaches 
a degenerate critical point in the direction of a centre manifold. 
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